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ABSTRACT 

Context. HH 158, the jet from the young star DG Tau, is one of the few sources of its type where jet knots have been detected at 
optical and X-ray wavelengths. 

Aims. To search, using Very Large Array observations of this source, radio knots and if detected, compare them with the optical and 
X-ray knots. To model the emission from the radio knots. 

Methods. We analyzed archive data and also obtained new Very Large Array observations of this source, as well as an optical image, 
to measure the present position of the knots. We also modeled the radio emission from the knots in terms of shocks in a jet with 
intrinsically time-dependent ejection velocities. 

Results. We detected radio knots in the 1996.98 and 2009.62 VLA data. These radio knots are, within error, coincident with optical 
knots. We also modeled satisfactorily the observed radio flux densities as shock features from a jet with intrinsic variability. All the 
observed radio, optical, and X-ray knot positions can be intepreted as four successive knots, ejected with a period of 4.80 years and 
traveling away from the source with a velocity of 198 km s"' in the plane of the sky. 

Conclusions. The radio and optical knots are spatially correlated and our model can explain the observed radio flux densities. 
However, the X-ray knots do not appear to have optical or radio counterparts and their nature remains poorly understood. 

Key words. ISM: jets and outflows - Herbig-Haro objects - radio-continuum: ISM - stars: individual: DG Tau 



1. Introduction 

] HH 158, the jet from DG Tauri was first reported by Mundt & 
Fried (1983), who presented an Ha image showing a well de- 
fined HH knot at ~ 8" to the SW of the star connected to DG Tau 
itself by a faint bridge. High resolution spectroscopy of this jet 

. was presented by Mundt et al. (1983) and Solf & Bohm (1993). 

' In the latter paper, the jet emission was traced to within ~ 0''2 
from DG Tau. 

DG Tau is located in the sky approximately in between the 
L1495 region and the star HP Tau. There are accurate distance 
determinations from very long baseline interferometry geomet- 
ric parallax to both the L1495 region (131.5 pc; Torres et al. 
2007; 2009) and HP Tau (161 pc; Torres 2009). Here we adopt 
for DG Tau a distance of 150 pc, intermediate to those of L1495 
and HP Tau. 

The proper motions of the knots observed up to distances of 
~ 10" were derived by Eisloffel & Mundt (1998), using several 
frames obtained over a ~ 7 yr time span, giving velocities of ~ 
150 km s"' (assuming a distance of 150 pc). Comparing adaptive 
optics images obtained with a ~ 2 yr time base, Dougados et al. 



(2000) obtained proper motion velocities of ~ 200 km s"' for 
the knots within ~ 6" from the source. These velocities imply a 
dynamical timescale of ~ 40 yr for HH 158 at that epoch. 

McGroarty & Ray (2004) suggested that two groups of HH 
knots (HH 702 to the SW and HH 830 to the NE) at angular 
distances of ~ 11' from DG Tau might be associated with the 
same outflow. However, McGroarty et al. (2007) showed that 
proper motion measurements only support the possibility of HH 
702 being associated with the outflow from DG Tau (the HH 
830 knots have proper motions which are not aligned with the 
outflow axis). The HH 702 knots have proper motions of ~ 240 
km s"', resulting in a dynamical timescale of ~ 2000 yr 

Spectroscopic data with 2D angular resolution (but with 
lower spectral resolution) of the region around DG Tau were pre- 
sented by Lavalley et al. (1997, 2000). These observations show 
that some of the ejections from DG Tau have a remarkable, bow 
shock-like morphology, and also show in a clear way the pres- 
ence of a faint counterjet to the NE (marginally seen in some of 
the previous observations). 
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Interestingly, the only optical spectrophotometric study of 
HH 158 with an extended wavelength coverage appears to be the 
one of Cohen & Fuller (1985). The spectrum described by these 
authors covers from ~4000 to ~7000A, and shows (reddening 
corrected) ratios of [O III] 5007/H/3=0.28, [O I] eSOO/Ho-^OJO, 
[N II] 6548+83/Ha=0.64 and [S II] 6716+31/Ha=0.69. These 
observed [O III]/Hy8 and [S ll]/Ha ratios identify HH 158 as a 
high excitation HH object (see Raga et al. 1996). 

More recent work has focussed on the low excitation lines. 
The paper of Solf & Bohm (1993, mentioned above) and the 
more recent papers of Bacciotti et al. (2000) and Coffey et al. 
(2007, 2008, who present red and near UV STIS spectra), and 
Pyo et al. (2003, who study the [Fe II] 1.644 jim emission) do 
not cover the blue region of the spectrum. Because of this, the 
[O III] 5007 emission reported by Cohen & Fuller (1985) has 
not been re-observed. 

Herczeg et al. (2006) obtained FUV spectra of DG Tauri, 
reporting the detection of fluorescent H2 lines, and the non- 
detection of lines like CIV 1549, which would be expected in 
the spectrum of a high excitation HH object. This non-detection 
possibly indicates that the high excitation emission region de- 
tected by Cohen & Fuller (1985) might be absent two decades 
later, or that it lies outside the region sampled by the slit in the 
STIS spectrum of Herczeg et al. (2006). 

The high excitation nature of HH 158 has been confirmed 
by the somewhat surprising detection of extended X-ray emis- 
sion along the outflow (Gudel et al. 2005, 2007, 2008, 2011; 
Schneider & Schmitt 2008; Gunther et al. 2009). These obser- 
vations show an X-ray knot at ~ 5" from DG Tau, along the 
direction of the HH 158 flow. 

We have obtained multiepoch VLA images and a red [S II] 
image of the region around DG Tauri, in order to explore the 
cuiTent morphology of the outflow, and to be able to relate the 
X-ray emission to the structures observed at other wavelengths. 

The base of the HH 158 flow was detected in the VLA obser- 
vations of Cohen & Bieging (1982) and Bieging et al. (1984). In 
the present paper we present new VLA observations made at 3.6 
cm in 1994, 1996, and 2009. We use these data, together with 
the knot positions measured over the past 20 years at optical, IR 
and X-ray wavelengths in order to derive a kinematical model 
for the evolution of the HH 158 outflow. 

The paper is organized as follows. In section 2, we describe 
the new and archival optical and radio continuum observations. 
In section 3, we present the sequence of radio continuum maps, 
and the red [S II] image. In section 4, we derive a simple kine- 
matic model for the time-evolution of the knot structure of HH 
158 over the past 20 years, and a model for the free-free con- 
tinuum produced in variable ejection velocity jets. Finally, the 
results are summarized in section 5. 



2. Radio and optical observations 

2.1. Radio continuum observations 

Since the objects in Taurus are known to exhibit relatively large 
proper motions (tens of mas yr"' , i. e. Loinard et al. 2007; Torres 
et al 2007) and our observations cover large time intervals (15 
years), we considered it necessary to first determine the proper 
motions of DG Tau to coiTect the positions and understand better 
any possible changes. For this proper motion determination we 
used the three epochs observed by us and described in detail be- 
low, as well as five additional observations taken from the VLA 
archive. Finally, we also included an EVLA (Expanded Very 



Large Array) observation taken in 2011 February 26 at an av- 
erage wavelength of 5.3 cm (two 1 GHz bandwidths centered at 
4.56 and 7.43 GHz) in the B configuration, as part of the Gould's 
Belt Distance Survey (Loinard et al. 201 1). The full set of 9 ob- 
servations covers about 30 years. The radio positions of DG Tau 
as a function of time are shown in Figure [T] The proper motions 
of the star derived from the fits shown in this Figure are: 

fia cosd = -1-7.5 + 0.9 mas yr"' 
Us = -19.0 + 0.9 mas yr '. 

These proper motions are consistent within 2-cr with those 
reported by Ducourant et al. (2005) from optical observations. 
Since our proper motion determination is about a factor of two 
more accurate than that of Ducourant et al. (2005), we adopt it 
and use for the position of DG Tau the following values: 

Qf(2000) = 04'' 27"' 04:'6880 + 0:'00055 x (epoch - 2000.0) 
5(2000) = +26° 06' 16"011 - 0"0190 x (epoch - 2000.0), 

where epoch is the epoch given in decimal years. In the images 
discussed below we use this position corrected for the epoch of 
the observation. 

The Very Large Array (VLA) observations made by us in 
the continuum at 3.6 cm were taken in three epochs spanning 15 
years. The parameters of these observations are summarized in 
Table 1 . The data were edited and calibrated using the software 
package Astronomical Image Processing System (AIPS) of the 
USA National Radio Astronomy Observatory (NRAO). 

In Figure|2]we show the images of the 1994 and 1996 epochs, 
both made in the highest angular resolution A configuration. The 
images are strikingly different. The 1994 image shows a single 
component slightly elongated in the NE-SW direction. The total 
flux density of this source is 0.41 + 0.04 mJy. In contrast, the 
1996 image shows two components, one is similar to that seen 
in 1994 but there is a new component 0!'42 + 0"01 to the SW 
of the first. If we assume that this new component is a knot that 
was ejected between the two epochs of observation and a dis- 
tance of 150 pc, we derive a lower limit of 1 16 + 3 km s"' for its 
velocity in the plane of the sky. The flux densities of these com- 
ponents are 0.44 + 0.03 mJy (NE component, coiTesponding to 
the star) and 0.40 + 0.03 mJy (SW component, corresponding to 
the knot). 

In addition to the high angular resolution observations made 
in 1994 and 1996, we recently made a deep integration of the 
region also at 3.6 cm but in the configuration C, that provides an 
angular resolution of ~ 3" to search for extended components 
around DG Tau. The image obtained from these observations 
is shown in Figure [3] Two components dominate this image: a 
bright one with a flux density of 1.47 + 0.02 mJy coincident 
with the star and a fainter one at 6"98 + 0"40 to the SW of the 
star and with a flux density of 0.15 + 0.02 mJy. We attribute this 
second source to a knot ejected in the past by DG Tau. 

2.2. Optical observations 

In order to see the present optical structure of the HH 158 out- 
flow, we have obtained an image of this object in the night of 
February 23, 2010. The narrow band image of DG Tau was ob- 
tained at the 2.6m Nordic Optical Telescope (NOT) of the Roque 
de Los Muchachos Observatory (La Palma, Spain) using the 
Service Time mode facility. The image was obtained with the 
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Frequency 


On-source 


Number of 


Phase 


Bootsttapped 


Beam 


Epoch/Configuration 


Project 


(GHz) 


Time (min) 


Antennas 


Calibrator 


Flux Density (Jy) 


Angular Size 


1994 Apr 16(1994.29)/A 


AR277 


8.46 


206 


18 


0403+260 


0.739±0.002 


0:'19x0:'16; -68° 


1996 Dec 24 (1996.98)/A 


AR277 


8.46 


631 


14 


0403+260 


0.691 ±0.002 


0:'25 X 0:'24; + 19° 


2009 Aug 13 (2009.62)/C 


AR694 


8.46 


414 


26 


0403+260 


2.34±0.02 


3:'07x2:'93; -15° 




EPOCH 




EPOCH 

Fig. 1. Position of the radio emission of DG Tau as a function of time 
from VLA and EVLA data. The right ascension (top) is with respect 
to 04'' 22™ and the declination (bottom) is with respect to +26° 06'. 
The dashed lines are least-squares fits to the data that give the proper 
motions discussed in the text. 



Andalucia Faint Object Spectrograph and Camera (ALFOSC) 
in imaging mode. The detector was an E2V 2Kx2K CCD with 
a pixel size of 13.5/im, providing a plate scale of 0.19 arcsec 
pixel-'. A [S II] filter (central wavelength A = 6725 A, FWHM 
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O 



04.75 04.70 04.65 04.60 

RIGHT ASCENSION (J2000) 



26 06 17.5 



17.0 



1996.98 




04 27 04.80 04.75 04.70 04.65 04.60 

RIGHT ASCENSION (J2000) 

Fig. 2. VLA contour images of the 3.6 cm emission associated with 
DG Tau for the epochs 1994.29 (top) and 1996.98 (bottom). The images 
were made with the weighting parameter ROBUST = 0. The contours 
are -3, 3, 4, 6, 8, 10, 12, 15, 20, and 25 times 19 and 15 fiJy beam~', 
the rms noises of the 1994.29 and 1996.98 images, respectively. The 
crosses mark the position of the star, determined as described in the text. 
The beams are shown in the bottom left comer and their dimensions are 
given in Table 1. 



= 60 A) was used to obtain an image of DG Tau in the [S II] 
6716, 6731 A emission lines. 

Two exposures of 900 seconds each were combined to obtain 
the final image. The angular resolution during the observations, 
as derived from the FWHM of stars in the field of view, was of 
0.9 -1.0 arcsec. The images were processed with the standard 
tasks of the IRAF reduction package. 
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Fig. 3. VLA contour images of the 3.6 cm emission associated with 
DG Tau for the epoch 2009.62. The image was made with the weighting 
parameter ROBUST = 5 to emphasize extended structure. The contours 
are -3, 3, 4, 6, 8, 10, 12, 15, 20, 30, 40, 60, 80, and 100 25 times 8.8 
yujy beam"', the rms noise of the image. The cross marks the position 
of the star, determined as described in the text. The beam is shown in 
the bottom left comer and its dimensions are given in Table 1. 



The [S II] image is shown in Fig. |4] In this figure we also 
show the positions found through paraboloidal fits to the posi- 
tions of the source and of the two observed knots. As the star 
is saturated, we have actually elliminated the central region, and 
carried out the fit to an unsaturated ring around the center of the 
PSF. 

Through these fits, we deduce distances of 6"75 and 12!'92 
from the source to the two knots along the jet. The same po- 
sitions are recovered in the two exposures to within 0!'01. 
However, while the first knot shows a well defined peak, the sec- 
ond knot (further away from DG Tau) shows a diffuse structure 
with a width of ~ 3". Because of this lack of a well defined 
peak, we do not use the second knot for calculating proper mo- 
tions (see section 3). 

3. Comparing optical, radio, and X-ray knots in DG 
Tau 

In Table 2, we give the positions of knots along the HH 158 jet 
compiled by Pyo et al. (2003), together with the positions of the 
radio continuum knots seen in Figs. 2 and 3, the position of the 
knot closest to the source in our optical images (see §2.2) and 
the X-ray knot positions of Gudel et al. (2008, 201 1). 

We have carried out a least squares fit to the knot positions x 
as a function of the time t of the observations of the form: 



X - xq + vjit -np) , 



(1) 



where the values of xq, vj and p are the same for all of the ob- 
served knots, and n is allowed to have values n = 0, 1, 2 and 3 
(choosing for the individual knots the values of n that minimize 
l\\Qx^)- This fit assumes ballistic (i.e. constant velocity) motions. 
From the fit, we obtain 



p = (4.80 + 0.32) yr; xq = (550 + 24)" 



vt = (0.277 + 0.012)"/yr = (198 + 9) km s" 




Fig. 4. Red [S II] image of the HH 158 outflow. The three crosses repre- 
sent the positions deduced from paraboloidal fits to DG Tau (upper left 
cross) and the two knots detected along HH 158. The image is displayed 
with a logarithmic scale (given in arbitrary units by the top bar). 



(2) 



where the velocity in km s"' was calculated assuming a distance 
of 150 pc to HH 158. This functional form is appropriate for a 
system of knots that travel with a constant velocity vj- (in " lyr 
for X in arcseconds and t in years), with an ejection period p. The 
values of n correspond to the successive knots. 

From the least squares fit to the observed knot positions, we 
divide the observed knot positions into four sequences (kO, kl, 
k2 and k3, corresponding to n = 0, 1, 2 and 3, see Table 2). 
Fig. |5] shows the four resulting linear x vs. t dependencies (cor- 
responding to n = 0, 1,2 and 3) together with the observed knot 
positions. 

It is clear that the observed knot positions can be interpreted 
as four successive knots, ejected with a period p - 4.80 yr and 
travelling away from the source at a velocity vj - 198 km s"' 
(see Eq. |2] and Fig. |5]l. This velocity is in good agreement 
with previously measured proper motions in HH 158 (see, e. g., 
Dougados et al. 2000). 

From Fig. 5 it is clear that the two radio knots reported by 
us are spatially correlated with optical knots that were observed 
close in time to the radio observations. On the other hand, the 
two observations of the X-ray knot (Giidel et al. 2008; 2011) 
show that it is the same knot that was observed as an optical 
knot close to the star in the early 1990's and has displaced sev- 
eral arc seconds since then. However, in the 2010.05 X-ray ob- 
servations the knot does not show optical or radio counterparts 
that should have beeb seen in the 2010.15 and 2009.62 obser- 
vations, respectively, presented by us. This lack of counterparts 
adds to the puzzling nature of the X-ray knot in HH 158, whose 
emission requires models with shock velocities between 400 and 
500 km s ' (Giinther et al. 2009) that do not appear to be present 
from data at other wavelengths. 
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Table 2. Optical, radio and X-ray knot positions 



distance ["] 



Obs. time 


kO 


kl 


k2 


k3 


ref. 


1991.24 




0.25 






1 


1992.86 


2.25 








2 


1994.84 


2.7 


1.4 


0.1 




3 


1996.98 






0.42 




ra 


1997.04 


3.3 




0.6 




4 


1998.07 


3.6 




0.93 




5 


1999.04 






1.3 




6 


2000.95 








0.45 


7 


2001.92 








0.75 


8 


2005.6 




4.32 






9 


2009.62 


6.98 








ra 


2010.05 




5.5 






10 


2010.15 


6.75 








op 



ra: radio, this paper 
op: optical, this paper 
1: optical, Kepner et al. (1993) 
2: optical, Solf & Bohm (1993) 
3: optical, Lavalley et al. (1997) 
4: optical, Dougados et al. (2000) 
5: optical, Lavalley et al. (2000) 
6: optical, Bacciotti et al. (2000) 
7: optical, Takami et al. (2002) 
8: optical, Pyo et al. (2003) 
9: X-ray, Gudel et al. (2008) 
10: X-ray, Gudel et al. (2011) 




1990 



2000 
t [yr] 



2010 



Fig. 5. Positions of the knots along the HH 158 jet as a function of 
time. The filled circles indicate the observed optical positions, while 
the open circles and the stars indicate the observed radio and X-ray 
knot positions, respectively. The four straight lines represent the least 
squares fit described in section 3 (the four lines corresponding to the 
values n = 0, 1,2 and 3, from left to right). 



4. A model for the radio-continuum emission from a 
bipolar outflow 

Raga et al. (1990) developed analytic and numerical models for 
jets from sources with intrinsically time-dependent ejection ve- 
locities. These authors showed that supersonic variabilities in the 
injection velocity in a supersonic flow result in the formation of 
two-shock structures (called working surfaces) that travel down 
the jet. More recently. Canto, Raga & D'Alessio (2000) devel- 
oped a method for solving the equations of a supersonic outflow 
with time-dependent parameters (ejection velocity and mass loss 
rate), based on considerations of momentum conservation for the 



internal working surfaces. In particular, these authors obtained 
solutions for a sinusoidal velocity variability with constant mass 
injection rate, and with constant injection density. In another 
work, Gonzalez & Canto (2002) developed a model to explain 
the observed free-free emission at radio frequencies around low 
mass stars. In their model, the ionization is produced by inter- 
nal shocks in a spherical wind themselves produced by periodic 
variations of the injection velocity. In this section, we present a 
model for calculating the radio-continuum (free-free) emission 
from a stellar flow with conical symmetry. We assume that the 
radiation is produced by internal working surfaces which move 
inside the bipolar outflow. 

4.1. Dynamics of a working surface 

We consider an outflow which is expelled with constant mass 
injection rate m, and with an injection velocity V(.(t) of the form. 



Vf (t) - V,, - Vc sin(wT), 



(3) 



where v„ is the mean velocity of the flow, is the amplitude 
of the velocity variation and o) is the angular frequency of the 
variation. 

From the formalism developed in Canto et al. (2000), it can 
be shown that the first working surface (in each cone) is formed 
at a distance from the source given by. 



rc = -- 
to 

at a time. 



v,v [1 - (vc/v„) sin(wTc)]^ 



(Vc/Vw) COS LiJTc 

^ _ ^ _ J_ [1 - (Vc/vw) sin(fa;Tc)] 
" w (vc/v„) cos(wTc) 

where r^. is the ejection time, 

1 . Vl + 8(Vc/v^)2 



n 
a> 



sin 



2(vv/v„) 



(4) 



(5) 



(6) 



Defining the variables f - (ti + T2)/2 and At = (t2 — Ti)/2, 
being ti and the ejection time of the flow directly downstream 
and upstream of the working surface, respectively, then the equa- 
tion that describes the dynamical evolution of the working sur- 
face is. 



a\r sin'^(6>f) + b\r sm(cof) + Cai- = 0, 
where. 



a^r = 



— KcoArf - sm\coAT)], 



b^r = {oj\t) [sin(wAT) - {(o\t) cos(coAt)], 



C\t = — sin (wAr) [sin (coAt) — (wAr) cos (coAt) 



(7) 



(8) 



(9) 



— {(oAt) sin (mAt)]. 



(10) 



The -I- sign in the solution of equation ^ gives the solution with 
|sin(6;T)| < 1. Using At as the free parameter (in the interval 
[0, ;r/£L>]), we obtain f by solving equation ^ and then we find 
Ti and T2. We then use the formalism presented in Canto et al 
(2000) to calculate the position and velocity of the working sur- 
face. 
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Fig. 6. Position r^j and velocity v^j of the working surface formed in 
a stellar jet with a sinusoidal variation in its velocity at injection. The 
parameters of the model are v„ = 300 km s"', Vc = 200 km s"', and 
<u = 1 .26 yr. We have assumed a constant mass loss rate of M = 5x10"* 
Mq yr~' . The physical description of the plots is given in the text. 



In Figure|6] we show an example for the position and veloc- 
ity of the working surface formed in a stellar jet with a sinusoidal 
ejection velocity variability. In our model, we have assumed a 
constant mass loss rate M = 5 X 10~* Mq yr~', mean velocity 
Vh, = 300 km s~', amplitude = 200 km s~*, and frequency 
w = 1 .26 yr~' (which corresponds to an oscillation period P ^ 
5 yr). It can be observed from the figure that the working surface 
is not formed in the flow instantaneously, but at a time tc = 2.67 
yr and at a distance = 26.81 AU from the central star Later, 
the working surface begins to be accelerated asymptotically ap- 
proaching the velocity v^. 

4.2. The geometric model 

We consider a stellar outflow expelled from a central star with a 
sinusoidal variation in the injection velocity and constant mass 
loss rate. We assume conical symmetry for the jet with an open- 
ing angle of the cones 0a, and an inclination angle between 
the outflow axis and the plane of the sky 0,-. Figure Q shows a 
schematic diagram of the bipolar ejection. 

A periodical variation in the ejection velocity (eq. ||3|) gen- 
erates a set of outgoing working surfaces. Assuming that the 
working surfaces do not lose mass sideways, we can apply the 
results presented in section 3. 1 to obtain their dynamical proper- 
ties (position and velocity as functions of time). It is possible to 
construct analytic (or semi-analytic) models of working surfaces 
in variable jets under one of the two following assumptions: 

- the working surfaces eject sideways all of the material pass- 
ing through the two shocks (associated with the working sur- 
face), so that the equation of motion is determined by a ram 
pressure balance condition, 

- most of the material going through the shocks stays within 
the working surface, so that a center of mass formalism can 
be constructed. 





' Bipolar outflow 




Mr axis 








■ / e. 

Workinp 




Surface 



To the observer 



Fig. 7. Schematic diagram showing a bipolar ejection with conical 
symmetry from a central star. The opening angle of the cones is 0^ and 
the inclination angle of the outflow axis [z'] with the plane of the sky 
ly, z] is 6i. The axes y = y' are perpendicular to the figure. 



If one compares predictions from (semi) analytic models (with 
these two assumptions) with axisymmetric numerical simula- 
tions, it is clear that the more realistic, axisymmetric solutions 
to the Euler equation lie in between the analytic models (see, 
e.g., Cabrit & Raga 2000). In general, the mass conservation as- 
sumption is prefeiTed because it leads to fully analytic solutions 
for many forms of the ejection variability (Canto et al. 2000), 
the ram pressure balance condition leads to an ordinary differ- 
ential equation that normally has to be integrated numerically 
(see Raga and Canto 1998). Actually, the models computed un- 
der the two approximations described above do not have stong 
differences, so that the choice of one over the other does not in- 
troduce important differences with more realistic jet models (i.e., 
axisymmetric numerical simulations). 

Under these assumptions, the first working surface forms at 
a time tc, the second one forms after one period P, and so on. In 
this case, the equation that describes the dynamical evolution of 
the working surfaces is. 



sin [{2k + 1) ;r - wr] = 



-bAr + (ft^^ - 4 «ArCAr) 
lUAr 



1/2 



(11) 



where ^ is a non-negative integer (0, 1, 2,...). The solution of 
equation ( fTTT i with ^ = gives the dynamics of the first working 
surface, and, in general, the solution with k = j — 1 gives the 
dynamical evolution of the j-th working surface. 

In order to obtain the flux density from the bipolar outflow, 
we must find the conditions that indicate us if a working surface 
is intersected or not by a given line of sight. In our simplified 
model, we assume that every working surface can be described 
as a portion of a sphere (a polar cap) whose physical size de- 
pends on the opening angle 0a and its position from the central 
source. Thus, let [x,y, z] and [x',y', z'] be two frames of refer- 
ence shown in Figure |7] Every point of the j-th working surface 
(at a distance r^jy] from the central star) satisfies the equation of 
a sphere in both reference systems, that is. 



ws[j\ 

and, therefore. 



. / 2 2 2^1/2 

X = ±{r ,.-3' -z) , 



(12) 



(13) 



where the symbol ± indicates the possibility that a line of sight 
intersects twice the working surface. 
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Therefore, the transformation equations between the two 
frames of reference are. 



x' = ±(r^^yj -y^- zH'^ cos 0,- - z sin 0,-, 



t! = ±(r^ r , - - z^f^ sin 0i + z cos 0,-. 

ws[j] 



(14) 



The intersection conditions of the j-th working surface, 
formed in the approaching (z' > 0) cone or the receding (z' < 0) 
cone, are obtained by comparison of equation [14] with the edges 
of the caps (±r^j|j] cos 9a)- As a consequence, the j-th working 
surface is intersected by a given line of sight when. 



z' > r^sU] cos Oa , 



or. 



z' < -r^sU] cos 0a. 



(15) 



(16) 



Assuming that the observer is located at a distance D from 
the source, a given line of sight intersects the plane of the sky 
[y, z] at the point {D sin sin O, D sin cos O), where and 
O are the inclination and the azimuthal angles, respectively. 
From equations (fT4l) -(fT6b. we obtain the intersection conditions 
in terms of these new variables. 



± (r^ , sin^0 sin^O - sin^0 cos^O)'^^ sin 6*; 



+ D sin0 cosO cos 0i > ry,s\j\ cos Qa > 
or, 

± , sin^0 sin^O - sin^0 cos^O)^^^ sin 



(17) 



optical depth of each working surface has the the contribution of 
the internal and external shocks. Using this representation, the 
optical depth of the j-th working surface is given by. 



Twsy] — 



y3, vllitj) no,i(tj)+pi v]'(tj) no,2itj) 



-2.1 



(19) 



where «o,i is the preshock density of the external shock, and 
no,2 is the preshock density of the internal shock at its time of 
dynamical evolution tj = t — (j — 1}P. 

At a given line of sight (specified by the angles 6 and 3>), we 
add the contribution of the i-th intersected working surfaces to 
obtain the total optical depth Tv(0, 0) along this line of sight, 
that is. 



Tv(0,O) = 



7=1 



(20) 



where fij = cos Oj, being Oj the angle between the line of sight 
and the normal vector to the j-th working surface at the intersec- 
tion point. It is easy to show that fij can be written as. 



1 + - — 0*1 - 1) 

wslf] 



1/2 



Finally, the flux density at radio frequencies from the bipolar 
outflow can be calculated by. 



(21) 



where By is the Planck function in the Rayleigh-Jeans approx- 
imation (By = 2kTe v^/c^ being k the Boltzmann constant, Tg 
the electron temperature and c the speed of light), and 0c = 
atan (r^sm/D). 



+ Z) sin0 cosO cos 0; < — r^sij] cos 0a , (18) 

where we have assumed that the observer is far enough that all 
points of the polar caps are located at the same distance from the 
observer (£) » r„s[j]). 

4.3. Predicted radio-continuum emission 

We consider the model described in section 3.2. First, we add 
the optical depths of the working surfaces intersected by each 
line of sight to obtain the total optical depth along this line of 
sight. Then, we use this optical depth to estimate the intensity 
emerging from this direction. Finally, the total flux emitted by 
the system can be estimated by integrating this intensity over the 
solid angle. 

Using the numerical models developed by Ghavamian & 
Hartigan (1998) for the free-free emission for a planar inter- 
stellar shocks, Gonzalez & Canto (2002) estimate the average 
optical depth of a shock wave. Assuming an average excitation 
temperature of lO^'K, these authors found that their results can 
be represented by Ty = j8no v~^"^, where no is the preshock 
density, Vj the shock velocity, and v is the frequency. The con- 
stants j8 and y depends on the shock speed. We note that the 



4.4. Numerical example 

In this section, we present a numerical example for the predicted 
radio-continuum flux at /l = 3.6 cm from a bipolar outflow with 
a sinusoidal ejection velocity. The opening angle of the cones 
is 6a = 30° and the inclination angle between the outflow axis 
and the sky plane is 0,- = 42°. The outflow is ejected from the 
central star with a mean velocity v„= 300 km s~^, an amplitude 
Vc = 200 km s~' and with an oscillation period P = 5 yi {a> = 
1.26yr~^). We have assumed a constant mass loss rate m = 5 X 
10~* M© yr~^ and a distance D = 150 pc from the observer 
These values are consistent with the parameters estimated from 
HH 158 and other jets from low-mass young stars (Eisloffel et 
al. 2000; Pyo et al. 2003; Agra-Amboage et al. 201 1). 

In Figure [8] we present our results. Once the first working 
surfaces (in both cones) are formed, at a time t = 2.67 yr (see, 
also. Fig. |6]l, the flux increases reaching a maximum value of 
"S'8.46 GHz ~ 0.5 mJy at t ^4 yr Afterwards, the emission de- 
creases until new working surfaces are formed in the outflow, 
and the flux increases again. We note that the flux density shows 
a periodical behaviour with the same period as that of the injec- 
tion velocity variability. Finally, we note from the figure that the 
predicted values of the flux density are in agreement with the 
radio observations of HH 158 reported in section 2. 
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Fig. 8. Predicted radio-continuum fluxes at /I = 3.6 cm for a bipolar 
outflow with a sinusoidal ejection velocity. The assumed parameters are 
v^= 300 km s~', Vc = 200 km s"', and oj = 1.26 yr"^. We have adopted 
a constant mass loss rate m = 5 x 10"* Mq yr~', and a distance D = 
150 pc from the observer. The opening angle of the cones are 0„ = 30° 
with an inclination angle S,- = 42°. The physical description of the plot 
is given in the text. 

5. Summary and Conclusions 

We presented an analysis of archive as well as new Very Large 
Array observations of DG Tau that detect emission from knots 
in the jet associated with this star. Radio knots are detected in 
the 1996.98 and 2009.62 data and they are found to correlate 
with optical knots in observations made close in time to the radio 
observations. One of these optical observations is provided by us 
in this paper. In contrast, the X-ray knot that is observed in the 
Giidel et al. (201 1) data does not coincide with the radio/optical 
knot and appears to be part of a different, later ejection that was 
first detected as a moving optical knot in the early 1990's. 

All the observed knot positions (optical, radio, and X-ray) 
can be interpreted as four successive knots, ejected with a period 
p = 4.80 yr and travelling away from the source at a velocity 
Vt = 198 km s~^. The next knot ejection is expected to take 
place around epoch 2014.0. 

We have modeled successfully the observed radio continuum 
emission in terms of working surfaces produced in a jet with a 
velocity at injection that varies sinusoidally with time. 
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